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1 Introduction 

The phenomenon of string breaking, as a consequence of confinement, is one of the most 
fundamental aspects of the gauge theories that lie at the basis of our understanding of 
high energy physics. While a comprehensive understanding in many cases, such as QCD, 
is still lacking, a lot of insight has been gained thanks to lattice gauge theory (LGT). 
In LGT, originally pioneered by Wilson [1], the theory is formulated on a discretized 
space-time lattice in such a way that the gauge symmetry is preserved and the continuum 
model is recovered in the limit of vanishing lattice spacing. This allows powerful Monte 
Carlo simulations which besides mass spectra [2] and phase diagrams [3] have been also 
successfully used to study static properties of confinement and of string breaking [4-9]. 
Despite the great success of this method, it suffers from the sign problem [10], which 
limits the kind of feasible simulations, and in particular real-time dynamics (although 
recently introduced techniques have allowed some real-time calculations in certain regimes 
[11, 12]). It is then of undeniable interest to explore other techniques that can overcome 
the limitations of standard LGT techniques. 
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A particularly promising approach is the application of tensor networks (TN) to the 
Hamiltonian formulation of LGT [13-24], These methods are free from the sign problem 
and allow the computation of real-time evolution [25-27], therefore opening up a new per¬ 
spective of addressing dynamics of problems that are currently intractable with Monte 
Carlo simulations. The first numerical results for static properties [14-16, 21], finite tem¬ 
perature calculations [22, 23] and out-of-equilibrium dynamics [16, 21] indeed show the 
power of the method. 

Another direction explored in recent years is quantum simulation of gauge theories [17, 
28-47]. The proposals for quantum simulators typically work as well with the Hamiltonian 
formulation of the LGT and map the gauge theory to the physical degrees of freedom 
of another quantum system which can be controlled experimentally thus being free from 
purely numerical limitations. Although a fully fledged quantum simulator is not available 
yet, this route seems promising for the future as the first proposals might come into reach 
with nowadays high level of experimental control. 

In this work, we consider a non-Abelian lattice gauge theory with dynamical charges 
in one spatial dimension. The Hamiltonian we simulate realizes an exact SU(2) gauge 
symmetry using a finite-dimensional representation for the bosonic degrees of freedom [47, 
48]. This model is hence a suitable candidate for the design of an atomic quantum simulator 
of SU(2) LGT. The finite dimension of the link Hilbert spaces allows the simulation of 
the model via direct application of MPS techniques, which we employ here in order to 
numerically explore the statical and dynamical aspects of the string breaking phenomenon. 

In order to detect and characterize the breaking of the string, we propose three different 
observables. They can be monitored during the MPS simulation, and should in principle 
also be accessible in an experiment, so that our results will be measurable in a potential 
future quantum simulator. 

More specifically we study three different scenarios. First, we consider the static as¬ 
pects of string breaking by determining the ground state of the system, which includes fully 
dynamical fermions, in the presence of two additional external static charges. This calcu¬ 
lation relates closely to the ones traditionally accessible by lattice Monte Carlo methods. 
We show that we can reliably determine if the string is present, and therefore identify the 
regions where we expect string breaking to occur. These calculations additionally demon¬ 
strate the suitability of the proposed observables for the detection of string breaking in 
dynamical scenarios. In the second place, we investigate the dynamics of string breaking 
by introducing the external static charges on top of the interacting vacuum and evolving 
the state in real time. When the string breaks, we can explicitly observe the screening of 
the charges via the creation of dynamical particles. Finally, we analyze how the picture 
changes when the charges added to the vacuum are themselves dynamical, a scenario which 
is closer to more realistic out-of-equilibrium situations. Also in this case we can recognize 
the string breaking if the fermion mass is small enough. 

The rest of the paper is organized as follows. In section 2 we review the model and 
describe the numerical methods we are applying, with a special focus on the MPS techniques 
we are using. In section 3 we present our numerical results for the static calculations. 
Section 4 contains the results for the real-time dynamics of string breaking with two static 
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external charges. In section 5 we present the results for the real-time dynamics of string 
breaking when a pair of dynamical charges is added to the interacting vacuum. To conclude, 
we summarize our findings in section 6. 


2 Model and Methods 


2.1 Truncated theory 

The model we consider is the Hamiltonian formulation of a 1+1 dimensional SU(2) LGT 
with dynamical fermions, in which the gauge symmetry is exactly realized with finite¬ 
dimensional link variables [47]. The model can be understood as a truncation of the com¬ 
plete SU(2) gauge theory, formulated as the Kogut-Susskind Hamiltonian with staggered 
fermions [49] 


H =£^ (''/4 t/ n''/++i + h-c.) 

n 

+ m ^(-l) n V4^n + 

n n 


( 2 . 1 ) 


where if> n = 



is a two component spinor taking into account the two “colors” of 


fermions (“red” and “green”) on site n. U n is a matrix on color space in the fundamental 
representation. Its entries, Un , act on the link between the fermionic sites n and n + 1 and 
change the corresponding electric flux. The operator is a group scalar, as it does not 
carry any color index, and gives the flux energy on link n. Physical states, ]<(>), have to be 
eigenstates of Gauss Law, i.e. G"|<^>) = (?"]</)) for all sites, n, where 


G" = L" - 72"^ - Q“, a E {x, y, z}. 

Q" = is the non-Abelian dynamical charge at site n. where cr" are the usual Pauli 

matrices. A nonzero value for q° indicates a non-vanishing static external charge. 1 * * L", 72" 
are the generators of left and right gauge transformations acting on link n, which correspond 
to the left and right electric held on the link and fulfill the commutation relations [L", Lm\ = 
[72", Rm\ = iSnm'E.y e “^ 7R n and [L",724] = 0. They are related to 
the operator for the electric flux energy as = Yla^n^n = Yla ^n^n- The Gauss 
Law components do not commute among themselves, as [G",G4] = — iS nm . e"^ 7 Gn, 
and thus cannot be diagonalized simultaneously. However, they all commute with the 
Hamiltonian, [77, G"] = 0, so that the Hilbert space is a direct sum of sectors characterized 
by their configuration of external charges, •{+"} [47]. 

In the strong coupling limit, corresponding to e = 0, the ground state is given by links 
carrying no flux, odd sites occupied with two fermions, and empty even sites. This state 

1 We are not considering these external charges explicitly in the Hamiltonian. As they are static, they 

would only contribute a constant mass term to eq. (2.1) leading to a constant offset in energy in all our 

simulations, thus being negligible. 
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can be written in a product basis as 

l^sc) = |2)|0)|0)|0)|2)|0).... 

In the previous expression bold numbers represent the occupation number of a fermionic 
site, while |0) represents the state of a link with no flux. The fermionic part of this 
state corresponds to the Dirac sea, and it is easy to check that it fulfills Gauss Law with 
G° |0sc) = 0 Vn, a. Applying on this state the gauge invariant string operator 

Snl = ^ n U n . . . Un+l-i^n+l (2.2) 

or its adjoint, S ^ nl , for odd l , 2 generates a particle-antiparticle pair at sites n and n + l, 
connected by a flux tube of length l. Such configurations will have an excess of energy of 
2m due to the particle-antiparticle pair, plus a flux energy proportional to the string length, 
l. Consequently, from a certain length on, it will be energetically favorable to reduce the 
flux energy by creating extra particle-antiparticle pairs, leading to configurations with a 
broken string. 

In the complete SU(2) model, the flux on a link is not bounded, and the dimension of 
the Hilbert space for each link is infinite. It is nevertheless possible to consider a theory 
truncated in a gauge invariant manner [48], where the maximum flux a link can carry is 
limited and the Hilbert spaces of the links are finite dimensional. Following the method 
in Ref. [48], each matrix element Un can be decomposed as a sum over all irreducible 
representations, which may be separated to summands that are gauge invariant themselves. 
Here we consider the model corresponding to the simplest non-trivial truncation of the full 
theory, meaning that only the trivial and the fundamental representation are kept, resulting 
in dimension 5 for the links. 

2.2 Numerical Methods 

We use the MPS ansatz to describe the system [50]. A MPS with open boundary conditions 
for N sites is given by 

l*>= E 

where the A l £ are complex matrices in <C DxD for 1 < k < N and A* 1 (A JJV ) is a row 
(column) vector. The states |*fc)^ =1 form a basis of the (^-dimensional local Hilbert space 
on site k. The number D, called bond dimension of the MPS, determines the number of 
variational parameters and limits the amount of entanglement that can be present in the 
MPS. 

The string breaking phenomenon can be studied statically or dynamically. Different 
well known numerical algorithms exist to find MPS descriptions of stationary or time- 
dependent states [26, 51]. In this work, we employ simulations of time evolution [25, 26, 52] 

2 As we are working with a staggered formulation the odd (even) sites correspond to antiparticles (par¬ 
ticles), consequently l has to be odd to create a string between an antiparticle and a particle. 
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where the evolution operator is approximated via a first order Taylor expansion. This allows 
us to study both static and dynamic scenarios by respectively using imaginary or real time. 
Additionally, the Taylor expansion preserves the symmetries of the Hamiltonian, such as 
gauge invariance. Then it is possible to explore a specific sector of the external charge 
distribution without explicitly implementing symmetries in the tensors. 

The numerical simulations have three main sources of error, each of them controllable 
by a suitable choice of parameters. The first one is due to the approximation of the 
evolution operator via a Taylor series. This error can be controlled by choosing the time 
step suitably small. Another source of error is due to the limited bond dimension. By 
obtaining results with different bond dimensions, the size of this truncation error can also 
be estimated and controlled. As we are working with finite systems, a third source of error 
arises from finite size effects which can be avoided by using sufficiently large systems. 

Although TN and in particular MPS can be formulated in terms of fermionic degrees 
of freedom, we choose to translate the fermionic degrees of freedom in Hamiltonian (2.1) 
to spins via a Jordan-Wigner transformation for our numerical simulations (details about 
the procedure and the relevant operators in spin formulation can be found in the appendix 
A) and in the following we work with the spin formulation. 

2.3 Detection of string breaking 

Detecting the presence or absence of the flux string during the evolution requires ob¬ 
servables that are suitable for the dynamical case. The Wilson loop and more general 
correlation functions, widely used in lattice calculations to determine the static potential, 
are typically evaluated in the limit of large Euclidean time and therefore not suitable for 
our setup [5, 8, 53]. Instead, to detect strings and string breaking in the dynamical setup, 
we propose three different observables. We monitor their values throughout the evolution, 
taking advantage of the fact that we have access to the MPS wave function at all times. 

In the first place, the spatially resolved spin and flux configurations in the system allow 
us to visualize the change with respect to the initial configuration (see e.g. figure 3). 

A second observable is the local imbalance between red and green species, also spatially 
resolved. A string operator (2.2) changes the fermionic content of the sites at its end-points, 
such that, on the strong coupling vacuum, it will produce a superposition of states having 
a single red or a single green fermion at the beginning or at the end of the string. This can 
be detected by the operator Q" 2 = \{n r ^ n — n g:7l ) 2 , where n nn and n g>n are the occupation 
numbers for the two fermionic species on site n. 3 

A third observable can be proposed that looks only at the flux content of the links. 
Applying a string of length l, starting at site n, on the strong coupling vacuum produces 
a state in which links between sites n and n +1 carry non-vanishing flux, whereas outside 
the region there is no flux. We can construct projectors P n i on this kind of configurations. 
However, we are not interested in a single string but rather in the statistics of string lengths 

3 One should note that the index a in Q“ 2 is not summed. This might look puzzling at first, as it seems 
that we are using a non color-neutral object as observable. However, a little calculation shows (see appendix 
A) that Q“ 2 is identical for all a. Therefore summing a would only yield an additional factor of 3 which 
we are not taking into account. 
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in a (time dependent) state | cf>(t)). Thus, we bin all strings of a certain length together 
and normalize by the number of possible strings of length l: 

' N-l 

In this way we can obtain histograms for the distribution of the string lengths in the state 
at a given time t. 

As we show in the next section, these three observables allow us to reliably determine 
whether our system still contains the initially imposed string or the string is broken. 

3 Ground state with static external charges 

The usual way of probing for string breaking in lattice Monte Carlo studies is the analysis 
of the static quark-antiquark potential. Using the MPS method, we can also determine 
the energy of an extra pair of external charges as a function of their distance by simply 
computing the ground state in a sector, where two external static charges are placed at the 
desired separation. 

In particular, we choose two static external charges q y = ±1/2, located in the central 
region of the system (to minimize finite-size effects) at a distance l. Gauge invariance 
requires that there is a flux tube connecting them. In the strong coupling vacuum the charge 
at each site is zero, as well as the left and right electric field on each link. Consequently, 
we can prepare a state with the desired external charge distribution by applying the (non¬ 
invariant) operator U n ... U n+ i_i on the strong coupling ground state. This operator only 
acts on links from n to n + l — 1, on which it creates a finite electric field. Hence the Gauss 
Law at the beginning and the end of the string yields a non-zero value corresponding to a 
state with qn = ±1/2 and q y n+l = =f 1/2. Subsequently we evolve this state in imaginary 
time to determine the ground state of (2.1) for a chosen set of parameters, (e, m, g). For the 
results presented in this section we use a time step At = 1.0 x 10~ 3 and a bond dimension 
D = 100, parameters which turn out to be sufficient to avoid noticeable numerical errors 
(see appendix B for a more detailed error analysis). 

In figure 1 we compare the ground state energy E for e = 3.0, g = 1.0 and various 
masses, where we subtracted the energy of the interacting vacuum without external charges 
£vac- For m = 3.0 and m = 5.0 we observe three regions as functions of the string length. 
For short strings the energy grows linearly with the string length, indicating the stretching 
flux tube between the charges, namely the presence of the string in the ground state. From 
a certain value l c on, the ground state energy does not depend on l. This is the signature 
for string breaking, as after reaching the threshold for creating particle-antiparticle pairs, 
reducing the flux is energetically favorable, the string breaks and the energy is independent 
of the initial string length. Finally we observe a third region, when the string length is 
already close to the system size, and finite size effects become noticeable. In our plots we 
clearly see that the values of l c , where one transitions from the string region to the breaking 
region, are independent of the system size. For m = 10.0 the mass is large enough that 
one does not leave the linear scaling region, even if we create the longest string that fits in 
our system. Hence, we do not expect string breaking to occur in this system. 
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Figure 1. Ground state energy for N = 22 (red crosses) and TV = 30 (blue X’s ) as a function of 
the initial separation of the external charges for m = 3.0 (left panel), m = 5.0 (central panel) and 
m = 10.0 (right panel). 

The difficulty of numerically detecting string breaking with Monte Carlo techniques has 
been traced back to the mixing of two-nreson states being hard to capture by Wilson loops 
[5, 6, 54]. With our method we can explicitly see this mixing happening, as illustrated by 
figure 2. The figure shows the energy as a function of imaginary time for several cases in the 
breaking region. At the beginning of the imaginary time evolution, we observe the energy 
going down until it reaches a metastable plateau, corresponding to the interacting string 
state, evidenced by the spin and flux configurations at that point, shown in the left inset 
panels. For later times there is then again a significant decrease in energy corresponding to 
the breaking of the string. This can also be seen in the spin and flux configuration (right 
inset panels of figure 2), where the region of high flux between the static charges is going 
away after the decrease in energy, and only a peak in the flux around the external charges 
remains. To better compare the string breaking scenario with that of a string ground 
state, the whole imaginary time evolution of spin, flux and charge square configurations 
for masses m = 3.0 and 10.0 and system sizes N = 22 and 30 is shown in figure 3. The flux 
configuration is giving an indication that the particle-antiparticle pairs created during the 
breaking process are clustering around the heavy external charges and screen the electric 
flux. Furthermore figure 3 reveals that there is essentially no difference between both 
system sizes. 

In figure 4 we plot the charge squared Q“ 2 and the histograms Pi for the initial con¬ 
figuration and for the final ground state respectively. We subtract the configuration of the 
interacting vacuum in the sector without external charges, in order to better visualize the 
difference. As one can see, in the breaking case two peaks around the external charges 
form in the charge square configuration, thus verifying that the particles created during 
the breaking process indeed cluster around the external charges and screen the flux. By 
contrast, in the nonbreaking case the charge square configuration only changes slightly. 
The histograms for the string lengths show a similar picture. In the nonbreaking case the 
clear initial peak at l = 11 is preserved, whereas for m = 3.0 it vanishes and peaks emerge 
around smaller string length. 

From figure 1 we can identify the parameter regions in which we expect string breaking 
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Figure 2. Energy as function of imaginary time for N = 30, m = 3.0, l = 11 (red solid line) and 
l = 21 (blue dashed line). The insets show the difference of expectation values for cr* n (red lines, 
left y- axes), a gn (green lines, left y- axes), and (blue lines, right y- axes) with respect to the 
interacting vacuum, shortly before and after the drop in energy (left upper panel: l = 21, t = 0.75, 
right upper panel: l = 21, t = 1.25; left lower panel: l = 11, t = 0.5, right lower panel: l = 11, 
t = 1.25). 

to occur. To study the real-time dynamics of the string breaking we select two distinct 
situations. We choose l = 11, which is deeply in the breaking region for m = 3.0 but still 
far enough from the point where finite size effects are noticeable, and for which we do not 
expect any string breaking with m = 10.0. For all the following real-time cases we set 
e = 3.0 and g = 1.0, as in the imaginary-time setup. 

4 Real-time evolution with static external charges 

MPS methods give us access not only to the static properties, but also to the real-time 
dynamics of the system. We may then investigate how the string breaking process manifests 
dynamically when the static charges are introduced in the interacting vacuum. To this 
end, we first compute a MPS approximation to the interacting vacuum of the theory using 
variational energy minimization [51, 55, 56]. Starting from this state, we apply again the 
non gauge invariant operator U n ... which effectively creates two static external 

charges q y = ±1/2 separated by a distance l = 11, in a gauge invariant manner, thus 
starting the connecting flux tube. Subsequently we evolve this state in real time with 
At = 1.0 x 1C± 4 , D = 100 and determine the time dependent spin, flux and charge square 
configurations along the chain, as shown in figure 5. As in the imaginary time case, one 
can see that a system size N = 22 is sufficient to avoid finite size effects. In order to better 
appreciate the dynamics, we show the details of the three proposed observables at fixed 
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Figure 3. Imaginary-time evolution of a string with length l = 11 between two external static 
charges, in a system of size N = 22 (upper three rows) and TV = 30 (lower three rows), for 
dynamical fermion masses m = 3.0 (left column) and 10.0 (right column). Shown are the site 
resolved expectation values for tr* , a gn (first and fourth row), J\ (second and fifth row) and Q“ 2 
(third and sixth row) as a function of imaginary time. 


times in the evolution, t = 0, 0.25, 0.5 and 2, for each mass in figures 6 and 7, where we 
again subtracted the interacting vacuum configuration to visualize the difference to the 
ground state without external charges. 

For m = 3.0 the plots in figure 6 reveal that between t = 0 and t = 0.25 particle- 
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Figure 4. Charge square configuration (upper panels) and histograms of string lengths (lower 
panels) at the beginning and the end of the evolution for a system of size N = 22 and Z = 11, the 
vertical lines in the upper panels indicate the position of the external charges. Panel (a) and (e) 
show the data for m = 3.0 at the beginning, panel (b) and (f) for m = 3.0 at the end, panel (c) and 
(g) for m = 10.0 at the beginning and panel (d) and (h) for m = 10.0 at the end of the evolution. 


antiparticle pairs are created inside the string region, leading to a significant increase in 
Q" 2 there. This is accompanied with a considerable decrease of the initial peak in the string 
length histogram around l = 11 and a change in the spin and flux configurations. This 
growth is continuing up to t = 0.5 where there has been a large amount of particles created 
in the string region and the peak in the histogram at l = 11 is already gone. At later 
times, t = 2.0, these particles are clustering around the region of the external charges and, 
as the flux configuration reveals, the external charges are screened, leading to a reduction 
of flux in the center region of the original string. For m = 10.0 the picture is significantly 
different, as there is essentially no change during the evolution in either the charge, spin 
and flux configurations or the string length histograms, which show a single dominant peak 
at l = 11 at all times, therefore indicating that the initial string is preserved during the 
evolution. The minor changes over time present in figure 7 result from the fact that the 
starting state of the evolution is not an eigenstate of the Hamiltonian, and consequently it 
is not perfectly steady. 

5 Real-time evolution with dynamical charges 

The scenarios with external static charges allow us to isolate and study the string breaking 
phenomenon. However, it is also worth investigating the alternative scenario in which the 
charges added to the vacuum are themselves dynamical, as this is arguably a more realistic 
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Figure 5. Real-time evolution of a string with length l = 11 between two external static charges, in 
a system of size N = 22 (upper three rows) and N = 30 (lower three rows), for dynamical fermion 
masses m = 3.0 (left column) and 10.0 (right column). Shown are the site resolved expectation 
values for er* ra , a* n (first and fourth row), J); (second and fifth row) and Q“ 2 (third and sixth row) 
as a function of time. 


situation. This enables additional configurations, in which the charges can move, to also 
play a role in the evolution, so that the string breaking phenomenon may be displayed 
differently. 

Using the same MPS techniques, we can also explore this setup. Thus we repeat the 
simulations described in section 4 but applying the gauge invariant string operator from eq. 
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Figure 6. Real-time snapshots for a system with external static charges at a distance Z = 11, 
dynamical fermion mass m = 3.0 and size N = 22. The upper row shows the charge square 
configuration, (Q“ 2 ) - (Q“ 2 ) v ac, the central row the spin, (a* /g n } - (cP /ffn ) V ac (red and green lines, 
left y- axes), and flux configuration, (J^) — (J 2 ) V ac (blue lines, right y-axes), and the lower row the 
histograms for the string lengths, P; — Pz jVac . The vertical red lines in the upper row indicate the 
position of the external charges. Each column corresponds to a time instant, t = 0, 0.25, 0.5 and 
2 . 0 . 

(2.2) on the interacting vacuum to construct our initial state. This again results in a state 
with a string between two charges. Different to the systems studied earlier, the charges are 
not external, but they are now created on a site and are fully dynamical. Again, we study 
the time evolution of the spin, flux and charge square configurations for different fermion 
masses and system sizes, as shown in figure 8, where we use At = 1.0 x 10 -4 and D = 100 
as in the previous section. We compare the case m = 10.0, in which the string does not 
break (see also figure 11), with m = 3.0. In the latter case, which exhibited clear breaking 
for static external charges, the fully dynamical situation shows differences (compare figures 
6 and 10). To better identify the features of string breaking, we look at an even smaller 
fermion mass, m = 1.0, shown in figure 9. 

In all cases we see a clear initial peak in the charge square configuration at the beginning 
and at the end of the string. For m = 1.0 and 3.0 new charges emerge especially in the 
string region and these peaks are quickly decaying. Also the flux in the string region is 
decaying while two small peaks are preserved roughly around the start and end point of 
the original string. The histograms for the distribution of string lengths show a similar 
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Figure 7. Real-time snapshots for a system with external static charges at a distance l = 11, 
dynamical fermion mass m = 10.0 and size N = 22. The upper row shows the charge square 
configuration, (Q“ 2 ) - (Q“ 2 ) v ac, the central row the spin, {<J*/ g>n ) - {&r/g,n)™c (red and green lines, 
left y- axes), and flux configuration, (J^) — (J 2 ) va c (blue lines, right y-axes), and the lower row the 
histograms for the string lengths, Pi — v ac- The vertical red lines in the upper row indicate the 
position of the external charges. Each column corresponds to a time instant, t = 0, 0.25, 0.5 and 
2 . 0 . 

picture: at the beginning of the evolution there is a clear peak at around l = 11 which is 
gone at around t = 0.25. In the m = 1.0 case the initial peak at l = 11 is less dominant 
as in the m = 3.0 case, due to the fact that in this case the interaction strength e is more 
important, leading to a state less close to a strong coupling string. For later times one sees 
that in both cases smaller string lengths are dominating in the system, therefore indicating 
that the string is broken. By contrast, in the case of m = 10.0 the dominant peak at l = 11 
is preserved during the entire evolution. Although the magnitude decreases, one can see 
from the flux and charge square configurations that there is not a lot of change, which 
indicates that the original string is still present. The slight changes have the same origin 
as in the previous cases. As we are not starting with an eigenstate of the Hamiltonian but 
rather do a local quench, the state is not perfectly steady. Furthermore in this case the 
charges are fully dynamical, what leads to richer dynamics. 
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Figure 8. Real-time evolution of a string with length l = 11 between fully dynamical fermions, in 
a system of size N = 22 (upper three rows) and N = 30 (lower three rows), for dynamical fermion 
masses m = 1.0 (left column), m = 3.0 (central column) and 10.0 (right column). Shown are the 
site resolved expectation values for <r* , cr* n (first and fourth row), (second and fifth row) and 
Q“ (third and sixth row) as a function of time. 


6 Conclusion 

We have studied the real-time dynamics of string breaking with MPS in a truncated 1+1 
dimensional SU(2) lattice gauge theory. By looking at the charge square, the flux config- 
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Figure 9. Real-time snapshots from the evolution of a string with length l = 11 between fully 
dynamical fermions, for m = 1.0 and system size N = 22. The upper row shows the charge square 
configuration, (Q% 2 ) - (Q“ 2 ) vac , the central row the spin, {<T*/ g>n ) - (K/ g ,n)^ (red and green lines, 
left y- axes), and flux configuration, (J^) — (Jj() vac (blue lines, right y-axes), and the lower row the 
histograms for the string lengths, Pi — P; !Vac . The vertical red lines in the upper row indicate the 
initial position of the charges added to the vacuum. Each column corresponds to a time instant, 
t = 0, 0.125, 0.5 and 1.5. 

uration and the statistics of string lengths, we are able to detect and characterize string 
breaking in statical and dynamical scenarios. 

We calculate the ground state in the sector with two static external charges. This 
allows us to clearly pick up the static signature of string breaking, and to identify parameter 
regions where we expect string breaking to occur. Furthermore our calculation shows that 
the string state is metastable during the imaginary time evolution, until configurations 
with a broken string mix in and reduce the energy, as expected. We have also shown that 
the particle-antiparticle pairs resulting from the breaking string cluster around the external 
charges and subsequently screen the electric field. This can be detected with the spin and 
flux configuration as well as in the charge square and the distribution of string lengths in 
the system. 

Using real-time simulations with MPS, we have studied the dynamics of string breaking 
in a setup with two static external charges. With the observables proposed before, the 
breaking and non-breaking cases are clearly distinguishable. We can explicitly observe 
that, in case the string is breaking, dynamical fermions are created that cluster around the 
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Figure 10. Real-time snapshots from the evolution of a string with length l = 11 between fully 
dynamical fermions, for m = 3.0 and system size N = 22. The upper row shows the charge square 
configuration, (Q“ 2 ) - (Q“ 2 ) v ac, the central row the spin, (<T*/ g>n ) - (^/ ffi „)vac (red and green lines, 
left y- axes), and flux configuration, (J^) — (J 2 ) V ac (blue lines, right y-axes), and the lower row the 
histograms for the string lengths, P; — Pz jVac . The vertical red lines in the upper row indicate the 
initial position of the charges added to the vacuum. Each column corresponds to a time instant, 
t = 0, 0.25, 0.5 and 2.0. 

external charges and screen them, therefore reducing the flux in the system. 

Finally, we have also simulated the time evolution of a string between fully dynamical 
fermions. We also identify situations in which the string breaks in this case. In particular, 
the decay of charges and the reduction of the flux in the middle of the string region indicate 
string breaking. Due to the fact that the initially created particle-antiparticle pair is now 
dynamical, we do not observe the clustering of the dynamically created charges, but rather 
a distribution along the system. 

In our study we have used a 1+1 dimensional SU(2) gauge model. This model could 
be a suitable candidate for the implementation of a quantum simulator for a SU(2) lattice 
gauge theory. The observables and setups that we describe would also be viable for such 
a potential experimental realization. Furthermore the applied techniques are not limited 
to the specific model studied here. The truncation method used from Ref. [48] works 
for arbitrary compact (or finite) Lie groups and is not limited to the case of one spatial 
dimension. Therefore also other gauge groups could be studied and with more general TN 
such as PEPS [57] this study could also be extended to higher dimensions. 
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Figure 11. Real-time snapshots from the evolution of a string with length l = 11 between fully 
dynamical fermions, for m = 10.0 and system size N = 22. The upper row shows the charge square 
configuration, (Q“ 2 ) - (Q“ 2 ) v ac, the central row the spin, - (^ /g n ) vac (red and green lines, 

left y- axes), and flux configuration, (J^) — (J^) vac (blue lines, right y-axes), and the lower row the 
histograms for the string lengths, P; — P; jV ac- The vertical red lines in the upper row indicate the 
initial position of the charges added to the vacuum. Each column corresponds to a time instant, 
t = 0, 0.25, 0.5 and 2.0. 


While completing the manuscript we became aware of a related work on TN and string 
breaking [58]. 
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A Spin formulation 

For the numerical implementation the fermionic degrees of freedom can be mapped to spins 
using a Jordan-Wigner transformation. Here we use the particular transformation from 
Ref. [59]: 
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The u-matrices are the usual Pauli matrices, where the subscript indicates on which color 
on site n they are acting. This labeling is convenient for practical purposes and completely 
equivalent to ordering the fermions as + 2 , ipg,2 ■ ■ ■ and doing a usual Jordan- 

Wigner transformation on a single fermionic field il> n which is related to original ones as 

Vh = VV,1, V>2 = Ipg, 1, + = +2, 

Using the transformation from eq. (A.l), the Hamiltonian can be expressed in spin 
language as follows: 
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The charge operators Q° = can then be expressed in spin language as 
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A simple calculation yields for the charge square components 
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a = x,y, z. 
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B Analysis of the numerical errors 


In this appendix we analyze the numerical errors of the results presented in sections 3-5. 
In the main text we have already presented data for different system sizes thus showing 
that N = 22 is enough to avoid noticeable finite size effects for the considered range of 
parameters. Here we focus on the influence of the time-step size and the bond dimension 
on our results. 

B.l Ground state with static external charges 

In figure 12 we show how the spin and flux configuration changes for a system of size 
N = 22, if a larger bond dimension or a smaller time step is used. As the figure shows, 
there is no significant improvement with larger bond dimension. Reducing the time step by 
a factor of two to At = 0.5 x 10~ 3 leads to a slightly changed spin and flux configuration on 
the order of 10 -2 for m = 3.0. However, these changes are predominantly at an early stage 
of the evolution and difference in the final configuration is much smaller. For m = 10.0 
the picture is qualitatively similar, at a very early stage there are small differences which 
are roughly one order of magnitude less than in the m = 3.0 case. All in all the changes 
for reduced time step and larger bond dimension are rather small compared to the data 
for At = 1.0 x 10~ 3 , D = 100 presented in figure 3 hence justifying our choice of time step 
and bond dimension. 

B.2 Real-time evolution with static external charges 

Also for the real-time evolution case with static external charges, we analyze our error 
analogous to the imaginary time case. The difference in results for a smaller time step 
At = 0.5 x 10 -4 and larger bond dimension D = 130 is presented in figure 13. A reduction 
of the time step shows a rather similar effect to the imaginary time case and leads to 
differences on the order of 10 -3 for m = 3.0 and 10~ 4 for m = 10.0. For m = 3.0, 
the differences comparing to results with enlarged bond dimension are bigger than in the 
imaginary time case, whereas for m = 10.0 there is essentially no change up to machine 
accuracy. Overall, we see the same picture as in the imaginary time cases that our choice 
of time step and bond dimension controls the error well enough to avoid considerable 
influences on the effects observed in figure 5. 

B.3 Real-time evolution with dynamical charges 

For the real-time evolution with dynamical charges the same error estimation as in the 
cases with static charges yields the results in figure 14. A reduction of the time step from 
At = 1.0 x 10~ 4 to 0.5 x 10~ 4 yields a change in the spin and flux configuration on the 
order of 10~ 2 . Contrary to the cases with static charges, the change is more pronounced 
for a large fermion mass m = 10.0, whereas for smaller fermion masses there is less change. 
For the bond dimension however, we see again the same effect that the changes are more 
pronounced if the mass is smaller, whereas there is almost no difference for larger mass. 
Nevertheless, also in this case the absolute changes in the spin and flux configuration are 
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Figure 12. Difference in the site resolved expectation values for spin (first and third row) and 
flux (second and fourth row) for a system of size N = 22 and dynamical fermion mass m = 3.0 
(left column) and 10.0 (right column). The upper two rows show the difference between results 
computed with a time step At = 1.0 x 10~ 3 and 0.5 x 10 -3 for D = 100, the lower two rows the 
difference between results computed with D = 100 and 130 for At = 1.0 x 10 -3 . 


rather small compared to the data presented in figure 8, thereby showing that the errors 
due to our choice of time step and bond dimension are negligible. 
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Figure 13. Difference in the site resolved expectation values for spin (first and third row) and 
flux (second and fourth row) for a system of size N = 22 and dynamical fermion mass m = 3.0 
(left column) and 10.0 (right column). The upper two rows show the difference between results 
computed with a time step At = 1.0 x 10 -4 and 0.5 x 10 -4 for D = 100, the lower two rows the 
difference between results computed with D = 100 and 130 for At = 1.0 x 10 -4 . 
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Figure 14. Difference of the site resolved expectation values for spin (first and third row) and flux 
(second and fourth row) for a system of size N = 22 and dynamical fermion mass m = 1.0 (left 
column), 3.0 (central column) and 10.0 (right column). The upper two rows show the difference 
between results computed with a time step At = 1.0 x 10 -4 and 0.5 x 10 -4 for D = 100, the lower 
two rows the difference between results computed with D = 100 and 130 for At = 1.0 x 10 -4 . 
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